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Abstract 


This study analyzes 16 molecular dynamic simulations of a biophysical model for B.-adrenergic (B2AR) and G protein-coupled receptor 
(GPCR) activation. In this model, a highly conserved cysteine residue, C106 (C3.25 or CysIII:01), provides a free sulfhydryl or thiol group in an 
acid-base equilibrium between uncharged (RSH) and charged (RS~) states that functions as an electrostatic molecular switch for receptor 
activation. The transition of C106 in the B2AR between acid and base states significantly changes the helical/transmembrane (TM) domain 
interactions and the electrostatic interaction energy differences (AAE"). The AAE™ changes correlate well with the experimentally observed 
ligand efficacies. The TM interaction energies display patterns compatible with those previously recognized as responsible for GPCR activation. 
Key differences between the agonist, epinephrine, and the antagonist, pindolol, are seen for the TM3 x 6, TM3 x 4, TM6 x 7 and TM1 x 7 
interaction energies. Pindolol also produces a weaker AAE"" interaction and less TM interaction energy changes, which are important differences 
between the agonist and antagonist ligands. The D115E mutant with pindolol displays a greater AAE™" and TM interactions than for the wild-type 
B2AR with pindolol. This explains the higher activity of pindolol in the D115E mutant. The constitutively active D130A mutant displays TM 
interaction patterns similar to those for the activating ligands implying a common pattern for receptor activation. These findings support the broad 
concept of protean agonism and demonstrate the potential for allosteric modulation. They also demonstrate that this two-state model agrees with 
many previous experimental and theoretical observations of GPCRs. 
© 2006 Elsevier Inc. All rights reserved. 
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1. Introduction 


Abbreviations: B2AR, f2-adrenergic receptor; DTT, dithiolthreitol; Receptor activation requires the recognition of an extra- 
D130A, the B2AR aspartate residue 130 to alanine mutation; D113E, the cellular signal that usually involves an endogenous agonist 
B2AR aspartate residue 113 to glutamate mutation; EPI, epinephrine; AEL, ligand activating its target receptor. The G protein-coupled 


the energy difference between the base and acid states for the electrostatic receptors (GPCRs), receive and transmit these signals to the 
energy terms; AAE™”, electrostatic interaction energy differences; GPCR, G 


protein-coupled receptor; G,, stimulatory GTP-binding regulatory protein of intracellular environment, As general models for all GPCRs, 
adenylylcyclase; ALJ, the energy difference between the base and acid states the B2-adrenergic (B2AR) and rhodopsin receptors have been 
for the Lennard—Jones energy terms; MD, molecular dynamics; NE, norepi- studied extensively to understand the complex molecular 


nephrine; PIN, pindolol; RMSD, root mean square deviation; ARH, the net changes that accompany receptor activation and signal 
change in a high affinity receptor state; TM, transmembrane 


Hs : : . . transduction. 
“ This work is dedicated to and in memory of Dr. Lester A. Rubenstein. : : : sete 
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B2AR receptors exhibit spontaneous activity similar to 
normally expressed receptors in the presence of an agonist 
ligand [1]. This observation separated receptor activation from 
the action of agonist ligands alone and prompted a revision of 
receptor models to include an intrinsically active receptor state. 
One consequence of this revision was that the resting 
populations of receptors must interconvert by themselves from 
resting to active states. However, the biophysical basis for these 
active and inactive receptor states has not been adequately 
defined or understood. One possible perspective is that there 
must exist an equilibrium between the inactive and active 
receptor states that can be perturbed to produce a shift in the net 
amounts of these states. An agonist ligand favoring the active 
receptor state would perturb the initial chemical equilibrium 
thereby inducing receptor activation in a manner similar to Le 
Chatelier’s principle. However, the underlying molecular 
mechanism of this perturbation needs to be more carefully 
examined. 

Mathematical models are essential for understanding the 
underlying biophysical events that determine the receptor 
response. In general, two-state mathematical models have been 
among the most successful for describing receptor activation 
[2-8]. Most of these models calculate either the proportional or 
fractional receptor occupancy as the overall receptor response 
[2-5]. Although it is seductive to assume that the proportional 
amount of an active receptor state should correlate with the 
biological response, the experimental evidence for receptor 
overexpression and spare receptors suggests that the calculation 
of the net change in the active receptor state is a much better 
measure for response than is the fractional or proportional 
change [6-8]. This is demonstrated by the effects of agonist/ 
antagonist combinations on the desensitization of B-receptors 
[7,8]. This is also demonstrated by receptors that are activated 
by overexpression since this requires a change between R and 
R” that is difficult to understand in terms of a proportional 
rather than a net change [7,8]. From this perspective, it is 
important to determine within the constructs of any molecular 
model what molecular states interconvert either by ligands, 
receptor overexpression or mutation. 

Several experimental techniques show that the relative 
movements of transmembrane domains II and VI (TM3 and 
TM6) correspond with receptor activation in both the rhodopsin 
and B2AR receptors [9-15]. Gether et al., using fluorescent 
labels, localized agonist-induced conformational changes of 
the B2AR to TM3 and TM6 [12]. Ballesteros et al. showed that 
charge-neutralizing mutations of E268 in TM6 and D130 in 
TM3 led to a significant increase in basal and pindolol- 
stimulated cAMP accumulation in COS-7 cells, which was 
presumably due to an increase in the relative motions of TM3 to 
TM6 [13,14]. In the rhodopsin receptor, Khorana’s group using 
spin labels attached to cysteine residues showed that rigid body 
movements of TM3 relative to TM6 were essential for 
activation [10]. 

Other investigators have also suggested that the cis—trans 
isomerization of retinal in rhodopsin causes repositioning 
followed by a movement of TM3 away from TM6 [9,11,15]. 
Clearly there are many studies that implicate TM3 and TM6 


with the molecular events leading to receptor activation. An 
implicit assumption from these important studies is that there 
exists an “off” conformation of the TM3 and TM6 domains 
that can be switched to an “on” conformation corresponding 
with GPCR activation. However, from this perspective, it is 
difficult to clearly translate TM3 and TM6 domain motions into 
discrete molecular entities for specific on and off states. 
Previously, we have attempted to bridge this gap by defining 
specific receptor states that relate the acid-base transition of a 
conserved sulfhydryl group to the mathematically defined R 
and R” states [16]. In this present work, we examine whether or 
not this model can be applied to the molecular dynamic analysis 
of the B2AR. 

Although controversial, there is a wealth of experimental 
evidence for the existence of at least one free sulfhydryl group 
that modulates the receptor response in many GPCRs [16-34]. 
For decades it has been known that there are free sulfhydryl 
groups in blood and on the external side of animal tissues in 
vivo [17]. However, most experimentalists and theoreticians 
have not seriously questioned that the conserved disulfide- 
bonded cysteines may be a potential source of free sulfhydryl 
groups under normal physiologic conditions. This may be 
primarily due to the fact that all of the GPCR families are 
considered to have a highly conserved disulfide bond 
connecting a conserved cysteine in the extracellular region 
of TM3 [19-21]. Within all of the GPCR families this disulfide- 
linked cysteine, C110 in rhodopsin or C106 in the B2AR, is the 
most conserved cysteine that is vital for GPCR expression and 
function [19-21]. This disulfide’s proximity to both the 
extracellular region and the ligand binding region suggests 
that the oxidation/reduction state of this bond could modulate 
both ligand binding and receptor function [23-25]. 

From specific sequence analyses, it is known that the B2AR 
receptor contains a total of 13 cysteine residues [18]. Of these 
cysteines, there are four potentially reactive cysteines, C106, 
C184, C190 and C191, on the external side of the B2AR that 
may form at least one extracellular disulfide bond [18]. 
Similarly, in the rhodopsin receptor there are three external 
cysteines, C110, C185 and C187 that are potentially reactive. 
The odd number of external cysteines in the rhodopsin receptor 
insures that there will be at least one free sulfhydryl group even 
when the disulfide bond is formed. 

The crystal structure for rhodopsin in the Protein Data Bank 
(1F88) shows C110 in the extracellular region of TM3 
positioned about 10 A away from the 11-cis-retinal chromo- 
phore and 7 A from glutamine E113 [15,22]. This C110 is 
normally considered to be in a disulfide bond with C187 and is 
apparently in a disulfide bond in the crystal structure; however, 
the nearby C185 apparently exists as a free sulfhydryl group, 
which may undergo acid-base transitions as well as other 
chemical reactions. Because the conditions for crystallization 
are very different from the natural state in which receptors 
normally function in vivo, it is difficult to imagine that the 
conditions for crystallization would preserve the free sulfhydryl 
groups of a reduced disulfide bond. Receptors transiently 
exposed to oxidizing conditions in vivo or in vitro may form this 
disulfide bond as a protective mechanism so that upon return to 
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reducing conditions, the normal receptor state could be easily 
regenerated. 

In this context, it is important to consider a few of the 
complex chemical equilibria between the acid/base and 
oxidized/reduced forms of cysteine [35]. First, the side chain 
can reversibly ionize: 


RSH > RS” +H* (1) 


Second, there is the oxidation/reduction equilibrium involving 
the disulfide bond: 


RSH + RS~ + RSSR + OH- (2) 


In the presence of a suitable oxidizing agent, a suitable base 
acts to increase the probability of disulfide bond formation by 
forming RS , which initiates the reaction cascade for disulfide 
bond formation (shown below) [36]. 

In general, many sulfhydryl groups or mercaptans are easily 
oxidized to disulfides and the interconversion between cysteine 
and cystine is important in biochemistry [36]. This mechanism 
has been studied for several oxidizing agents and varies with the 
agent. For oxygen it is: 


RSH + B` e RS +BH 


RST +0, = RS’ +0, 


RST +0: & RS’ +0,7 
2RS° + RSSR 


O27 +BH = OH“ +B +0, 


This mechanism with respect to the sulfur involves the loss 
of a proton, reduction to a free radical followed by radical 
coupling [36]. The resultant equation from this reaction cascade 
is Eq. (2) above. This reaction cascade demonstrates the 
important role played by base ions (B`) to favor free radical 
formation in the presence of an oxidizing agent such as 
molecular oxygen. This also demonstrates some of the 
experimental difficulties in studying free thiol groups since 
exposure to air in the presence of a base is often sufficient to 
oxidize free thiols to disulfides or several other oxidized species 
too numerous to mention here [35-39]. 

From the above considerations, assumptions that the disulfide 
cysteines are always resistant to chemical derivatization may not 
be correct for the normal conditions of receptor function in vivo. 
When one or more of the extracellular cysteines of the BZAR 
(C106, C184, C190 or C191) are either mutated or chemically 
altered thiol/disulfide exchange may occur with the remaining 
free cysteines or other free thiol groups to complicate the 
analysis. In addition, double mutants often fail to express either 
as an extracellular or functional form of the receptor, which 
makes it difficult to determine which of the four cysteines in the 
B2AR are most critical for receptor function and expression 
[18,19]. This is usually interpreted as evidence for the disulfide 
bond being essential for receptor activation [18,19]. Although 
this view has been accepted without very much critical 
questioning, there is at least one alternative explanation [16]. 


Evidence from experimental binding studies from several 
GPCRs demonstrate that a free sulfhydryl group is next to the 
agonist ligand binding site [26-33]. Sidhu made use of this fact 
to purify the D-1 dopamine receptor from rat striatum [26]. 
Others using a morphine ligand with an attached sulfhydryl 
group found a free sulfhydryl group near the ligand binding site 
of the u opioid receptor [27]. 

Additional studies have found free sulfhydryl groups next to 
the agonist binding site in the u, 6 and k opioid receptors 
[28,29]. Heinonen et al., using a series of SH-reactive 
“molecular yardsticks”, found a cysteine thiol group about 
8.6-14.5 A from the position of Asp113 (D113) in the az- 
adrenergic receptor [30]. Since this aspartate residue is 
important for attracting the positive amine group of binding 
ligands in many GPCRs including the B2AR, the proximity of 
free cysteines may be considered as a critical component of 
normal receptor function, which has been experimentally 
observed [23,24,31—34] and theoretically justified [16]. 

Experimental measurements of receptor activation suggest 
that the disulfide bond requires reduction to a free sulfhydryl for 
the functional state of several GPCRs in vivo [23—25,34]. One 
of the simplest and most direct experimental studies demon- 
strated that externally applied sulfhydryl reducing agents 
potentiate the response of tissues to various agonist ligands 
[34]. This study found that the reduction of disulfide bonds to 
free sulfhydryls increased the tissue responsiveness for many 
receptors and increased the potency of many agonists [34]. One 
problem with these studies is that physiologically opposing 
receptors may be simultaneously stimulated by the application 
of thiol reducing agents, which would confuse the analysis of 
the tissue response. 

Additional evidence for the requirement of reduction for 
normal receptor function is evident in a number of other 
important studies. A free sulfhydryl group in the active 
rhodopsin receptor state is suggested from the MI to MII 
equilibrium shift of 4 units from a pK of 4.1 for ligand-free 
opsin to 8.5 [3]. This is also supported by the findings of 
Saranak and Foster that implicate a free sulfhydryl in the 
reductive activation of the rhodopsin photoreceptor [23]. For 
the muscarinic receptors, Florio and Sternweis found that the 
reducing agent DTT produced functional activation and 
potentiated agonist activation [25]. Pedersen and Ross using 
purified B-adrenergic receptors and purified G, protein in 
phospholipid vesicles found that the treatment of these 
receptors with the sulfhydryl reducing agent dithiolthreitol 
(DTT) activated the G protein in either the presence or absence 
of agonist [24]. Since there was no effect of DTT on G, alone, 
they reasoned that it was the receptor itself that was the target of 
the reducing agent and not the G protein [24]. Several other 
studies have also shown that constitutive activation occurred in 
reconstituted systems containing both B-adrenergic receptors 
and G proteins when the membrane vesicles were prepared with 
sulfhydryl reducing agents [31-33]. These experimental 
findings need to be included into a comprehensive model for 
receptor activation. 

Previously, the acid and base states of the sulfhydryl group 
within the C148 residue in the 5-HT2A receptor were linked to 
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the inactive and active receptor states R and R“ [16]. In the 
current study, we have extended this model to the B2AR where 
the corresponding cysteine residue is C106 (C3.25 or 
CysIII:01). This cysteine, located on the extracellular side of 
TM3, functions as a molecular switch that toggles between the 
charged (RS) and uncharged (RSH) states associated with the 
acid to base transition of the sulfhydryl group [16]. This 
provides a specific molecular switch for receptor activation that 
correlates well with both the ligand induced electrostatic 
interaction energy difference (AAZ) and the experimentally 
measured ligand efficacies and the preferential attraction of 
ligands for the base receptor state (RS) [16]. In this context, 
ligand-induced perturbations produce a net change in the 
amount of the base receptor state (RS), that can mathema- 
tically describe the receptor response in a two-state model [6,7]. 

For this study, we examined 16 molecular dynamic (MD) 
simulations of full agonists, partial agonists and an antagonist 
of the wild-type and two mutant B2ARs. The D130A mutant 
changed the aspartate residue D130 in the DRY portion of TM3 
to alanine. This mutation removes a negatively charged residue 
from the B2AR and produces a constitutively active receptor in 
at least two separate systems [14,40]. The other mutant is the 
D113E that replaces the aspartate residue D113 (D3.32) with a 
glutamate residue. This mutation does not alter the net charge 
of the receptor, but profoundly alters ligand efficacies causing 
both the antagonist pindolol and the full agonist norepinephrine 
to behave like partial agonists [19]. This behavior is examined 
through these simulations of B2AR activation within the 
context of this two-state acid—base model for GPCR activation. 


2. Methodology 
2.1. The molecular model 


The B2AR model displayed in Fig. 1 was developed and 
further refined from the ‘‘What-IF” structures produced by 
Vriend at the EMBL labs in Heidelberg, which is based upon 
Baldwin et al.’s general model [20]. The D130A and D113E 
mutant structures were also modified from this B2AR model. 

For this model, the change of the sulfhydryl group in C106 is 
associated with the acid and base states in the B2AR (see 
Fig. 1), which allows the sulfhydryl group to exist in either acid 
or base states (SH or S~) [16]. Since the conserved Cys106 is 
located near the extracellular region, nearby water molecules 
would be expected to influence or modulate this electrostatic 
switch, but we did not model this possibility. 

The receptor models comprising the seven transmembrane 
regions were constructed with the program QUANTA. The loop 
regions were not considered in this particular model although 
future work may require the inclusion of these elements. All 
hydrogen atoms were included explicitly. Standard CHARMM 
charges were applied to all receptor atoms except for the base 
state of the cysteine (C106 or C3.25) in which case the charges 
were derived from a calculation of the molecule ethyl sulfide 
with a 6-31G basis set using the Hartree-Fock method. The 
calculated charges were then scaled and adapted to CHARMM 
format. Ligand models in the form of cations were built with the 


Extracellular 


C106 (3.25) 
em 


B2AR 


Intracellular 
(A) 


The B2AR Sequences: 


TM1(EC to IC)§: 
MGIVMSLIVLAIVFGN*VLVITAIA 


TM2(IC to EC): 
VTNYFITSLACAD*LVMGLAVVPFG 


TM3(EC to IC): 
106 113 130 
WCtGFWTSIDVLCVTASIETLCVIAVDR*YFAITS 


TMA(IC to EC): 
KARVIILMVWIVSGLTSFLPIQM 


TMS(EC to IC): 
IASSIVSFYVP*LVIMVFVYSRVFQEAKR 
TMG(IC to EC): 

KEHKALKTLGIIMGTFTLCWLP*FFIVNI 


TM7(EC to IC): 
ILLNWIGYVNSGFNP*LIYCRSPD 


B2AR mutants of the TM3 residues: 


D113E: (D3.32 or D113 to E113) 
D130A: (D3.49 or D130 to A130) 


§ (EC to IC) is the direction of this transmembrane domain from extracellular 
(EC) to intracellular (IC). Because of the way sequences are written, the direction 
of the EC and IC ends alternate. Therefore TM1, TM3, TMS and TM7 are EC to 
IC whereas TM2, TM4 and TM6 are IC to EC. 

* Conserved residues in each TM domain 1.50, 2.50, 3.50, 4.50, 5.50, 6.50 and 


7.50 (Ballesteros JA, Weinstein H, Methods Neurosci 1995; 25: 366-428). 


t C (C3.25 or C106 - in either the acid SH or base S state) 
(B) 


Fig. 1. (A) The B2AR model as viewed through the plane of the cell membrane 
(TM3 is the center solid red ribbon). A includes the initial docking position of 
epinephrine (EPI) as a representative ligand. Ligands were manually docked in 
the acid receptor state with no initial refinement of the side chains. The space- 
filled residue at the top is the highly conserved cysteine residue C106 (C3.25 or 
CyslIII:01), which provides a free sulfhydryl] or thiol group (yellow) in either an 
acid or base state (SH or S~). The space-filled residues at the bottom are the 
critical DRY region near the intracellular end of the receptor. (B) The sequences 
for this BZAR model are refined from ‘‘What-IF” structures proposed by Vriend 
(http://www.gpcr.org/7tm/models/vriend2/index.html) and based upon Baldwin 
et al.’s model [20]. The C106, in which the sulfhydryl group was modeled in 
either the acid, SH; or base, S`, state, is labeled in TM3. The two mutants, 
D113E and D130A, which occur in TM3, are listed below the sequences. 
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program QUANTA. Each ligand model was constructed with 
the chirality corresponding to the most active form. The models 
for epinephrine and norepinephrine had chirality R and the 
model for pindolol had chirality S. Charges were determined 
from calculations with the program GAUSSIAN-98 using 6- 
31G basis sets with the Hartree-Fock method. Energy 
minimization and molecular dynamics calculations were done 
by using the set of programs CHARMM (version 24). 

The ligand was docked into the acid form of the receptor 
with the cysteine in the acid state (SH) using facilities of the 
program QUANTA. The ligand was positioned near enough to 
form a salt bridge with either the aspartate D1 13 in the wild type 
receptor or the glutamate E113 in the D113E mutant receptor. 
The position of the ligand in the acid state was also used as the 
starting position for docking in the base state (S`). 

Subsequent calculations and molecular dynamic simulations 
were done with the C106 in either the uncharged SH (acid) or 
charged S (base) state. Initial receptor minimization included 
sequential intervals of steepest descent, conjugate gradient, and 
adopted basis Newton—Raphson. The dynamics calculations 
included heating to 310 K (30 ps), equilibration (50 ps), and 
simulation (200 ps) portions. The lengths of the bonds 
involving hydrogen atoms were constrained according to the 
SHAKE algorithm, allowing an integration time step of 
0.001 ps. Integration of Newton’s equation of motion was done 
by using the Verlet algorithm. 

Harmonic constraints with a force constant of 0.5 were 
applied in order to maintain the structure of the receptor 
backbone. A dielectric constant of 4.0 was used to emphasize 
electrostatic interactions between the ionized cysteine and the 
rest of the receptor since the electrostatic interactions computed 
generally take place within the protein interior in the absence of 
intervening solvent. A cut-off distance of 60 A was used to 
model the distant interactions between the ionized cysteine and 
the rest of the system. Equilibration was observed after 50 ps 
for all runs. For each run, the minimized average structures used 
for the analyses were taken from the last 5 ps of the simulations. 

Sixteen molecular dynamics (MD) simulations were 
analyzed for this study comprising eight acid—base pairs with 
the C106 cysteine in either the SH or S~ states. The MD 
simulations included the wild type (WT), and two mutant 
B2ARs with and without bound ligands as described in Table 1. 
These include each of the following receptor systems in acid 
and base states: the wild-type (WT) B2AR alone and with the 
ligands epinephrine (EPI), norepinephrine (NE) and pindolol 
(PIN); the D113E mutant alone and with the ligands NE and 
PIN; the D130A mutant alone. 

We understand that this B2AR model without loops, in 
vacuum, with backbone restraints may not completely validate 
a GPCR activation hypothesis, but hope that it will be 
considered in the larger context of the experimental and 
theoretical data that it attempts to explain. 


2.2. The calculation of interaction energies 


For each of the base and acid states of a particular receptor 
system the interaction energy (AE) was calculated for each 


Table 1 
CHARMM energy differences for the base—acid receptor states 
Receptor Ligand AE ALJ? AEL? AIF? 
(kcal/mol) 

System 
WT WT None +33.9 +10.4 +15.7 7.8 
WT + EPI WT EPI —3.8 -8.7 —5.1 10.0 
WT + NE WT NE —10.6 -9.0 —3.6 2.0 
WT + PIN WT PIN 47.7 +10.4 436 —6.3 

Mutants 
D113E D113E None +10.4 —5.9 +09 0.5 
D113E+NE D113E NE —8.3 —10.6 +0.6 1.7 
D113E+PIN D113E PIN +12.3 $5.0 +68 154 
D130A D130A None +10.6 +13.1 -O1 -2.4 


è ALJ =the energy difference between the base and acid states for the 
Lennard—Jones energy terms (kcal/mol). 

> AEL =the energy difference between the base and acid states for the 
electrostatic energy terms (kcal/mol). 

© ATE = the energy difference between the base and acid states for the internal 
energy terms (bending, stretching and torsional) (kcal/mol). 


molecular complex by taking the energy of the receptor—ligand 
molecular complex (Eg) and subtracting the energies of the 
free receptor and ligand molecules (Ep) and (EL) respectively: 


AE = Egy — Eg — EL (3) 


In general, the lower or more negative the interaction energy the 
more favored the formation of the molecular complex. 

Each of the interaction energies can be further partitioned 
into Lennard—Jones (LJ) and electrostatic (EL) terms according 
to the expression: 


AE = AEU + ABEL (4) 


The differences between interaction energies (AAF) for 
various molecular complexes were calculated as before [16] 
where the AE"" differences for the base minus the acid states 
(B — A) of a ligand-receptor complex is calculated as: 


AAE™ = AE™ (B) — AE™ (A) (5) 


All of the interaction energies were calculated from the 
average structures generated from the last 5 ps of the MD 
simulations. 


3. Results 
3.1. Differences between the acid and base states 


Fig. 2A and B displays the electrostatic potentials for the 
acid and base states of the B2AR plotted for the positive (blue) 
and negative (red) potentials. These figures all have the same 
orientation as the B2AR model displayed in Fig. 1. Fig. 2A and 
B show that the transition of C106 from the uncharged acid 
state (RSH) to the charged base state (RS) produces a change 
from an electrophilic (blue) to a more nucleophilic (red) state. 
This is consistent with the base state (RS ) preferentially 
attracting the positively charged cationic amine group of 
binding ligands and thereby creating a higher affinity state 
consistent with previous results [16]. 
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(B) 


Fig. 2. (A) Plot of the acid state with the same orientation as shown for the 
B2AR model in Fig. 1A with the extracellular end at the top of the picture and 
TM3 as the solid red ribbon. The space-filled residue at the top is the highly 
conserved cysteine residue C106 (C3.25 or CysIII:01), which provides a free 
thiol group (yellow) in an acid—base equilibrium between uncharged (RSH) and 
charged (RS) states. The space-filled residues at the bottom are the DRY 
region near the intracellular region of the B2AR. The electrostatic potentials are 
plotted as positive 25 (blue) and negative —25 kJ mol™' (red) meshes in both A 
and B. The acid state corresponds to the uncharged RSH or “R” receptor state in 
this model. (B) Plot of the base state with the same orientation as shown for the 
B2AR model in Fig. 1A and the acid state in A. The base state corresponds to the 
charged RS~ or “R”” receptor state. Portions of these composite pictures were 
generated with the Swiss-PdbViewer (from Ref. [51]). 


After performing the MD simulations representing a 
composite of several molecular structures that exist around 
the equilibrium energy, we analyzed various properties of the 
average structures. Fig. 3 displays the root mean square 
deviations (RMSDs) for the acid and base differences of the 
conserved residues in the wild type B2AR (labeled “*”) 


together with those residues displaying a RMSD greater than 
1A. Except for the arginine R3.50 (R131), the other conserved 
residues show less than 1 A difference. This suggests that these 
residues may function as a scaffold that is not as dynamically 
mobile as other parts of the B2AR. 

The distances between the alpha carbons of R131 and D130 
in TM3 and the alpha carbons of the residues within the other 
helices were generally shorter in the acid than in the base state 
suggesting that the base structure is somewhat larger than the 
acid structure (not shown). 

The distance between the R131 and D130 alpha carbons in 
the DRY region may also be affected by the strength of the salt 
bridge and the loss of this stabilizing interaction in the base 
state as discussed below. 

If the extra negative charge in the base state shows a higher 
affinity for the positively charged agonist ligands, then we must 
nevertheless explain why antagonist ligands do not favor the 
base over the acid receptor state. In order to understand this, we 
analyzed the acid-base interaction energy differences of the 
average molecular structures in more detail. 

The CHARMM energy differences for the WT, D113E and 
D130A B2AR mutants are displayed in Table 1 for the base 
minus acid difference (B — A). Apart from the extra negative 
charge on the thiol group, the base state of the WT B2AR has a 
higher overall CHARMM energy than the acid state (see 
AE(B — A), Table 1). 

In both the D113E and D130A mutant receptors, the base 
state also has a higher CHARMM energy than the acid state 
although the total energy change, AE(B — A), is less than that 
for the WT B2AR. The higher energy of the base state in the 
WT and the mutant B2ARs may be partly attributed to the 
repulsion energy between the negative charge on the ionized 
sulfhydryl group (RS ) and the other negatively charged 
regions such as the backbone oxygens and the carboxy] group 
of D113, which is two helical turns (about 10.8 A) away from 
the cysteine C106. Overall, the total CHARMM energy profile 
suggests that the acid state is energetically more favored than 
the base state in the unliganded native conformation, which is 
expected if the native receptor exists primarily in an off 
conformational state. 

For the agonist-bound receptors in almost all cases, the 
difference in the CHARMM energies favors the base state over 
the acid state, whereas the antagonist pindolol PIN showed an 
energy pattern that is contrary to this. This may be partially due 
to a larger ALJ(B — A) energy difference for PIN compared to 
the agonist ligands (see Table 1). 

For the D113E and D130A mutant receptors, there is a lower 
energy difference compared to the that for the WT receptor 
alone (compare AE(B — A) 33.9 with 10.4 and 10.6 kcal/mol in 
Table 1). Interestingly, the overall CHARMM energy differ- 
ences for these mutations are very similar but result from 
differing intramolecular forces. For the D130A mutant, the 
CHARMM energy of the AEL(B — A) difference decreases 
dramatically from 15.7 to —0.1 kcal/mol. In contrast, the D113E 
mutant shows CHARMM energy differences for the base minus 
acid states (B — A) that result from each of the separate energy 
terms ALJ, AEL, AIE and AE as shown in Table 1. These 
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Select TM Residues 


Fig. 3. The differences in the root men square deviation (RMSD) for the base minus acid states (B — A). The residues are listed as for the Ballesteros and Weinstein 
notation [52]. Conserved residues 1.50, 2.50, 3.50, 4.50, 5.50, 6.50 and 7.50 are starred “‘*”’. The other residues are either residues of interest such as the DRY region 


or those residues with the RMSD(B — A) > 1 A. 


observations suggest that there may be a lower energy barrier for 
the overall acid to base transition for these mutant B2ARs 
resulting from different molecular interactions. 

For the unbound wild-type receptor, the transmembrane 
(TM) domain interaction energies for the base minus acid states 
(B — A) are plotted in Fig. 4 for all possible TM pairs. These 
energy differences range from a low of —3 kcal/mol for 
TM2 x 3 to a high of 18 kcal/mol for TM6 x 7 followed 
closely by 16 kcal/mol for TM3 x 4. 

The largest total energy changes occur in TM6 x 7 and 
TM3 x 4 followed by TM3 x 6, TM1 x 7 and TM3 x 5 
(Fig. 4). The largest changes in the electrostatic energy (AEL) 
occur for the TM3 x 4 and TM6 x 7 interactions; whereas, the 
largest changes in the Lennard—Jones energy (ALJ) occur for 
the TM3 x 6 and TM6 x 7 interactions (Fig. 4). 

This suggests that these major helical interactions are 
perturbed by the transition from an acid to base receptor state 
apart from any ligand binding. TM7 also appears to be strongly 
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Transmembrane Pairs 


Fig. 4. Transmembrane domain interaction energy differences for the WT 
B2AR receptor displayed for the differences between the base and acid (B — A) 
states. The total energy differences for all TM pairs is displayed as the base 
minus acid (B — A) and is also broken down into the separate van der Waals 
(ALJ) and electrostatic (AEL) interactions. 


connected indirectly through the electrostatic interactions 
between TM6 and TM7. These interactions demonstrate the 
complex dynamics associated with the acid to base transition. 


3.2. Differences in interaction energies and their 
correlation with ligand efficacies 


The smaller ligands have a larger positive charge density and 
a more favorable interaction with the negatively charged 
cysteine in the base state suggesting that their relative efficacies 
may be related to their charge densities and their ability to 
interact with the base receptor state [16]. 

The CHARMM energy differences, AE(B — A), for the 
ligands EPI, NE, and PIN bound to the WT receptor are —3.8, 
—10.6, and +7.7 kcal/mol, respectively (Table 1). When 
norepinephrine and pindolol are bound with the D113E mutant 
receptor, the values of the CHARMM energy difference are 
—8.3 and 12.3 kcal/mol, respectively (see the AE@B — A) 
values for D113E + PIN and D113E + NE in Table 1). A partial 
explanation may be that for the NE and PIN ligands complexed 
with the D113E mutant receptor, the average distances between 
the ligands and the receptor increase in comparison with those 
in the wild-type complexes, which may be due to the slightly 
longer sidechain on the glutamate residue in the D113E mutant. 

The RMSDs were calculated for the acid and base states of 
the unbound receptor state and the ligand-bound receptor. The 
greatest RMSDs occur within TM6 and are greater for the full 
agonists, EPI and NE, than for the antagonist, PIN (not shown). 

However, we could not consistently connect the RMSDs of 
any specific residue to the experimental ligand efficacies. 
Although TM7 showed a larger RMSD with PIN in the D113E 
mutant than in the WT receptor (not shown), this may be 
partially attributed to the displacement of PIN by the larger 
glutamate residue, which places the ligand closer to TM7. 

In contrast, the electrostatic interaction energies (AAEH 
show that PIN in the D113E mutant receptor has a larger AAZ”™ 
than for PIN bound to the WT B2AR. One reason for this may 
be that the basic nitrogen of PIN is closer to C106, which 
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Table 2 
Electrostatic interaction energy differences for the base minus acid states with 
the percent sodium fluoride activity 


Table 3A 
Changes in key transmembrane domain interaction energy differences for the 
base minus acid states 


Receptor system AAE®" (kcal/mol) NaF” (%) 
WT + EPI —14.0 95 
WT + NE —14.4 95 
D113E+ NE —5.5 50 
D113E + PIN —6.4 35 
WT + PIN —4.7 0 


è These values were taken from Strader et al. [19]. The correlation coefficient 
for AAE®" vs. %NaF is r = 0.92. 


generates a larger AAE* in the D113E mutant. As a result, the 
electrostatic interaction energy of PIN with the negatively 
charged C106 residue is greater, which provides an explanation 
for the increased efficacy of pindolol as primarily due to the 
increased electrostatic favoring of the base state in the D113E 
mutant. 

The calculated AAE™ energies for WT+PIN were 
—4,7 kcal/mol versus —6.4 kcal/mol for the D113E+ PIN, 
which are in good agreement with the experimentally observed 
increase in the experimentally measured percent sodium floride 
activity as shown in Table 2. 

For the wild-type receptor with the NE ligand, the calculated 
AAE®" was —14.4 kcal/mol; whereas for the D113E mutant 
receptor with NE, the AAE*" was —5.5 kcal/mol suggesting a 
diminution in the efficacy of NE in the mutant consistent with 
previous experimental findings [19] (see Table 2). These 
findings are consistent with the observations from Strader’s 
group that the D113E mutation changes PIN from an antagonist 
to a partial agonist and reduces the activity of NE to that of a 
partial agonist [19]. 

In order to further understand the effects of ligand 
occupancy on the B2AR, we examined the transmembrane 
(TM) interaction energy differences for the acid and base states 
in the EPI bound WT receptor. In general, the binding of EPI 
changes the interaction energy differences of he base minus 
acid states (B — A) from positive to negative for the TM1 x 7, 
TM3 x 4, TM3 x 6 and TM6 x 7 interactions (see Fig. 5A). 
This suggests that one of the actions produced by EPI binding is 
to increase the magnitude of the TM interactions between 
TMI x 7,TM3 x 4,TM3 x 6and TM6 x 7 in the base state of 
the receptor. Many of these interactions have been previously 
linked either directly or indirectly to receptor activation [9-15]. 

With the WT B2AR bound EPI, the TM3 x 4 interaction 
energy changes from a positive 16 to a negative 23 kcal/mol 
suggesting that there is a closer interaction of TM3 and TM4 in 
the base state. The TM3 x 4 interaction is also the largest 
interaction measured for all of the EPI induced interaction 
energy changes shown in Table 3A. 

The interdomain interaction energy differences for the 
antagonist PIN in the WT receptor show an opposite pattern 
from the EPI and other activating ligands (compare Fig. 5A and 
F and Table 3A). Although PIN binding reduces the positive 
interaction energy differences for TMI x 7, TM3 x 4, 
TM3 x 6 and TM6 x 7, PIN fails to make them negative. 
This suggests that the weaker electrostatic interaction energy 


Receptor TM pair energies (kcal/mol) 
systems (B — A) 
TMI x7 TM3 x 4 TM3 x 6 TM6 x 7 

WT + EPI —7.3 —23 —7.7 —5.7 

WT + NE —3.2 —2.7 3.8 —6.1 
D113E+ NE 1.3 —0.5 2,2 —4.5 
D130A —2.7 —6.4 4.8 1.6 
D113E + PIN 1.7 —0.4 —1.3 1.7 
WT + PIN 1.4 T 3.8 5.4 


change (AAE*") for the PIN binding translates into TM 
interaction energy changes that fail to reach a critical threshold 
value for receptor activation. This may be a key difference 
between the antagonist and agonist ligands. 

For the mutant D113E + PIN receptor system, TM3 x 4 and 
TM3 x 6 become negative suggesting an activating pattern that 
is supported by previous experimental observations showing 
partial agonist activity [19,40]. From this analysis, TM3 x 4 
changing from positive to negative emerges as a consistent 
predictor for all of the activating ligands and mutations in this 
study. 

To further isolate the effects of ligand binding, we subtracted 
the TM energy differences of the WT(B — A) from the other 
receptor systems in order to see the net effect of the bound 
ligands. These results are displayed in Table 3B. This highlights 
those differences created by ligands or mutations that change 
the TM energy patterns of the unbound WT B2AR. 

By this analysis, the largest effect of EPI binding is more 
readily seen for the TM3 x 4, TM6 x 7, TM3 x6 and 
TMI x 7 interdomain interaction energies (Table 3B). These 
changes are largely consistent with the TM interaction energies 
in Table 3A, but more clearly demonstrate that the net effect of 
ligand binding centers around the major TM3 x 4, TM6 x 7, 
TM3 x 6 and TM1 x 7 interactions. 

In the receptor systems with the largest agonist activity, the 
TM3 x 4 and TM6 x 7 interactions are also the largest 
(Table 3B). The pattern of TM interaction energies for the 
mutant D130A receptor compared to the EPI bound WT 
receptor also displays a strong correlation (r = 0.90 in Fig. 6) 
suggesting that this may be a general pattern for both ligand 
and constitutive receptor activation. These general patterns 
support Christopoulos and Kenakin’s concept of protean 


Table 3B 
Changes from the wild-type B2AR in key transmembrane domain interaction 
energy differences for the base minus acid states 


Receptor TM pair energies (kcal/mol) 
ti B-A 

Sems Bm TMIx7 TM3x4 TM3x6 TM6x7 
EPI-WT —12 —39 —16 —24 
NE-WT —8 —19 —4 —24 
D113E + NE-WT —4 -17 —6 —22 
D130A-WT —8 —23 —4 —16 
D113E + PIN-WT -3 -17 —10 —16 
PIN-WT —4 —9 —4 —13 
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agonism by demonstrating that there exist similar activation 
patterns that vary around a common theme for activating 
ligands and mutations [41]. 

These findings also demonstrate that although pindolol is an 
antagonist, it contains some elements, albeit weaker, within its 
interaction pattern in common with the full agonist epinephrine. 
Parts of these patterns may represent energies that fail to reach a 
threshold energy change for receptor activation since the PIN 
energy differences in Table 3B are generally smaller than those 
for the agonist ligands. 

However, in the D113E mutant receptor, the antagonist 
pindolol changes to a partial agonist and norepinephrine 
activity decreases approximately 50% to that of a partial 
agonist. This occurs in our analysis by the alteration of key TM 
domain interaction energies such as the changes in the 
TM3 x 4 and TM3 x 6 for PIN (compare Fig. 5D and F) 
and those for the TM3 x 4 and TM6 x 7 for NE (compare 
Fig. 5C and B). 

In contrast, the accompanying changes to the TM domain 
interaction energy differences for D113E + PIN compared to 
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WT+PIN show a more negative energy pattern for the 
TM3 x 4, TM3 x 6 and TM6 x 7 interactions in the D113E 
mutant consistent with a higher activity for PIN in this mutant 
(Tables 3A and 3B). 

If the constitutively active mutant D130A is assumed to have 
approximately the same activity as the D113E+ PIN system, 
then the TM6 x 7 interdomain interaction energy differences 
listed in Table 3B correlate well with the experimentally 
measured relative intrinsic efficacies listed in Table 2 (% NaF). 
Interestingly, there is also a good correlation between all of the 
other ligand bound and mutant receptors with their experi- 
mentally reported activities (r=0.94 with the D130A, or 
r=0.92 without the D130A mutant). 


3.3. Regions within the receptor that may have a role in 
signal recognition and transmission 


In our model, C106 is two-helical turns (~10.8 A) above the 
important aspartate residue D113 that forms a counterion with 
the protonated amine of catecholamine ligands. This is 
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Fig. 5. (A-G) The transmembrane domain (TM) interdomain interaction energy differences for the base minus the acid (B — A) states for all possible TM domain 
pairs. In all parts, the vertical scales are equal. Each figure shows the patterns of the transmembrane domain (TM) interdomain interaction energy differences for the 
base minus the acid (B — A) states with the unliganded wild-type B2AR WT (B — A) in the black bars given for comparison. A shows epinephrine with the wild-type 
receptor (WT + EPI (B — A)). B shows norepinephrine with the wild-type receptor (WT + NE (B — A)). C shows norepinephrine with the D113E mutant 
(D113E + NE (B — A)). D shows pindolol with the D113E mutant (D113 + PIN (B — A)). E shows the D130A mutant (D130A (B — A)). F shows pindolol with the 
wild-type receptor (WT + PIN (B — A)). G shows the D113E mutant (D113E (B — A)). 
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consistent with Heinonen et al.’s, suggestion that a reactive 
cysteine or thiol group is 8.6-14.5 A from the Asp113 in the az- 
adrenergic receptor [30]. Similarly in rhodopsin, the distance 
from retinal to the glutamate E134 in the ERY region is between 
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Fig. 6. Correlation of the TM interdomain interacton energy differences of the 
base minus acid states for the constitutively active mutant D130A vs. the wild- 
type B2AR bound with epinephrine (D130A-WT(B — A) vs. EPI-WT(B — A); 
r = 0.90). 


19 and 30 Å (PDB 1F88) [22], which generally agrees with the 
distances for ligands docked within our B2AR model. 

The protonation and movement of R131 (R3.50) has been 
previously suggested as a potential mechanism for receptor 
activation [13,40,42]. Previously in the two-state acid-base 
model for the 5-HT2A receptor, the charged cysteine was 
observed to have long range electrostatic interactions with the 
aspartate and arginine residues within the DRY region [16]. 

In the present study, the acid to base change of C106 is 
accompanied by a corresponding displacement of the aspartate 
and arginine in the DRY region of the B2AR (0.7 and 1.7 A, 
respectively) (Fig. 3). This demonstrates that the acid to base 
transition perturbs R131 enough to disrupt this bonding. This 
suggests that deprotonation of C106 followed by a separate 
protonation of R131 may be a feasible molecular mechanism 
for the transmission of an extracellular signal to the 
intracellular side of the receptor. 

All of the possible pairs of TM domain interaction energy 
differences demonstrate complex patterns for both the 
electrostatic (AEL) and Lennard-Jones (ALJ) interactions as 
shown in Fig. 4 for the wild-type B2AR. The largest 
perturbations occur for the AEL energies of the TM3 x 4 
and TM6 x 7 interactions and the ALJ energies of the TM3 x 6 
and TM6 x 7 interactions. 
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For all of the MD simulations, the most consistent change 
accompanying receptor activation is the change of the TM3 x 4 
interaction energy from positive to negative together with the 
AAE™ becoming more negative. This is observed for all of the 
active systems we studied; WT + EPI, WT + NE, D113E + NE, 
D130A and D113E+ PIN, shown in Fig. 5A-E. These TM 
domain interactions may in turn affect the second and third 
intracellular loops ICL2 and ICL3), which are considered to be 
important in receptor-G protein-coupling [15,30]. These 
interactions may produce a repositioning of the arginine in 
the polar pocket surrounding the DRY region demonstrating 
that our relatively simple model generates complex electro- 
static and conformational changes within the B2AR that link 
TM changes to previously recognized patterns of receptor 
activation [5,13,21,40]. 


4. Discussion 


This work tests a two-state acid—base model for an ionizable 
cysteine residue found within all GPCRs. This cysteine is 
strategically located two-helical turns or about 10.8 A above 
the important aspartate residue D113 in TM3 that forms a 
counterion with the protonated amine of catecholamine ligands 
in adrenergic receptors [30]. The acid—base transition of this 
cysteine, C106 in the B2AR, modulates the internal energies of 
the receptor and the interaction energies of various ligands that 
interact with the B2AR in ways that correlate well with 
previously observed transmembrane domain interactions and 
the experimentally measured ligand efficacies (Table 2 and 
Fig. 6). Although it would be preferable to have a larger 
sampling of ligands with their experimental efficacies, these 
simulations taken together with previous studies suggest that 
those ligands that activate the B2AR also show a preferential 
electrostatic attraction for the base form of the receptor state 
[16]. 

Previous two-state models for GPCR activation have been 
considered to represent a simplification of the available 
conformational space [2,43]. However, single-molecule spec- 
troscopy of the B2AR receptor has recently suggested that there 
are two distinct states with a population enrichment of one state 
occurring when bound with agonist molecules [43]. The authors 
of this study suggested that there is a population of active 
receptor states that are separate from the active states in the 
native receptor, which appears contrary to the findings from 
Bond et al. [1]. 

One possible way to resolve these seemingly contradictory 
findings is through the interpretation that our two-state model 
provides that allows for a net enrichment of the active state 
through either receptor overexpression or ligand binding [6-8]. 

On the molecular level, we have determined that the discrete 
molecular states RSH and RS connect the mathematically 
general R and R” states to electrostatic interaction energy 
changes (AAE*") that closely correlate with their experimen- 
tally measured ligand efficacies as seen previously for the 5- 
HT2A receptor [16]. This suggests a molecular basis for the R 
to R” shift in the receptor equilibrium as a measure for the net 
receptor response and supports the concept that cysteine 


modulation of ligand binding may be an important mechanism 
to explain how ligands activate these receptors [16]. 

Although this picture may be more complex than previously 
characterized, the acid—base transition of C106 clearly creates 
energy perturbations that are largely compatible with many 
previous experimental findings [9—15]. In the unliganded wild- 
type B2AR, the TM1 x 7, TM3 x 4, TM3 x 5, TM3 x 6 and 
TM6 x 7 interactions are all strongly positive for the unbound 
WT (B — A) states of the B2AR (gray bars in Fig. 4 and the 
black bars in Fig. 5A—G). However, these interactions switch 
from positive to negative upon the binding of epinephrine with a 
rank order of TM3 x 4 > TM6 x 7 > TM3 x 6 > TMI x7 
(Fig. 5A). 

In addition, the TM3 x 4 interaction also changes from 
positive in the unbound D113E mutant to negative in the PIN 
bound D113E mutant, which suggests a possible threshold 
activation energy for this interaction (compare Fig. 5G with D). 
All of the MD simulations of the receptor systems experimen- 
tally shown to be active, WT + EPI, WT + NE, D113E + NE, 
D130A and D113E + PIN, consistently show that the TM3 x 4 
interaction energy changes from a positive to negative value 
(Table 3A). Although more work needs to be done to support 
this finding, at least one study has implicated TM4 in the 
activation of the M1 muscarinic receptor [44]. 

One of the more subtle changes involves the D113E 
mutation that changes pindolol from an antagonist in the WT 
B2AR to a partial agonist in the D113E mutant [19]. To our 
knowledge there have been no other theoretical models that 
have explained this experimental finding. However from our 
simulations, the change from an aspartate to glutamate creates a 
larger difference in pindolol’s interactions between the acid and 
base states, which is directly reflected in the AAE™ energy 
difference and the changes in the TM interaction energies. This 
is most evident by the TM3 x 4 and TM3 x 6 interactions, 
which change from positive for the WT + PIN to negative for 
the D113E + PIN system (compare Table 2 and TM3 x 4 and 
TM3 x 6 in Fig. 5D and F). 

Unlike the agonist ligands in the wild-type B2AR, pindolol 
fails to make the TM3 x 4, TM6 x 7, TM3 x 6 and TMI x 7 
interaction energies negative, which represents a crucial 
difference between the agonist and antagonist ligands. This 
molecular mechanism elucidates these subtle changes within 
the dynamic patterns of B2AR activation. 

Further experimental observations suggest that some ligand- 
dependent conformations favor one particular G protein over 
another as a basis for protean agonism [2]. Although EPI and 
NE are about equally effective at activating the B2AR, their TM 
interactions show some unique differences. For the NE bound 
B2AR, all of the major TM interactions except for TM3 x 6 
change from positive to negative; whereas for the EPI bound 
receptor the TM3 x 6 also changes from positive to negative 
(Fig. 5A and B). However, their overall TM changes show 
similar patterns and display a strong correlation relative to the 
acid-base changes of the WT receptor (r = 0.83). 

These common patterns are comprised of changes in the 
TMI x 7,TM3 x 4,TM3 x 6and TM6 x 7 interaction energy 
differences (Tables 3 and 4). Interestingly, the constitutively 
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active mutant D130A also displays a similar pattern of TM 
interactions compared to the EPI bound receptor suggesting a 
common pattern for receptor activation (Figs. 5A and E and 6). 
This supports the concept that protean agonism arises naturally 
from ligand or mutation specific patterns due to variations in 
internal interaction energies that may favor receptor binding to 
one particular G protein over another. 

In comparison with rhodopsin, Sakmar et al. implicated the 
H-bond networks linking the TM1 x 7, TM3 x 4, TM3 x 6 
and TM6 x 7 interactions to receptor activation [15] and 
Ballesteros et al. showed that the charge-neutralizing mutation 
of D130 in TM3 of the B2AR to asparagine (D130N) led to a 
significant increase in basal and pindolol-stimulated cAMP 
accumulation, which was thought to be due to an increase in the 
relative motions of TM3 to TM6 [13,14]. In the context of our 
model, these TM interaction changes arise from the acid to base 
transition alone. 

The “protonation hypothesis” for B2AR activation suggests 
that the protonation of the aspartic/glutamatic acid residue in the 
D/ERY portion of the cytoplasmic side of TM3 is a key event for 
initiating receptor activation [21]. The constitutively active 
D130A mutation is thought to free the arginine R131 or R3.50 
within the critical DRY region [14,40], which suggests that the 
protonation of the D130 or D3.49 frees this arginine to move out 
of the polar pocket [13]. This is also supported by the 
experimental finding that in isolated BZAR membrane prepara- 
tions decreasing pH favors receptor activation in vitro [45]. 

In the context of our model, the existence of at least two 
different pH-sensitive residues on both interior and external 
sides of the receptor may complicate these interpretations since 
another study has shown that an increase in external pH 
increases receptor activation [46]. This still raises the important 
question of how activating ligands that bind to the extracellular 
region might induce the protonation of acidic residues on the 
cytoplasmic side of the receptor. In our model, the aspartate and 
arginine within the DRY region move by 0.7 and 1.7 A (Fig. 3), 
which may be enough to disrupt an intramolecular bond such as 
a salt bridge. This would reconcile the observations that either 
an acidic or basic pH may promote receptor activation, since the 
differing pH-sensitivities of the two residues Cys 106 or Asp 
130 may be involved on both the extracellular and intracellular 
side of the receptor in the cell membrane. 

Recent experimental observations of receptor dimerization 
together with parallel receptor activation studies demonstrate 
that receptor dimer formation enhances receptor activation 
[47,48]. Zeng and Wess produced the first direct evidence for 
the existence of muscarinic receptor dimers formed from two 
conserved cysteine residues, C140 and C220 located on the 
extracellular side of the receptor with C140 corresponding to 
C106 of the B2AR [47]. 

In the dopamine D2 receptor (D2R), Javitch and co-workers 
showed that the site of crosslinking occurs at C168 on the 
extracellular end of TM4 and concluded that cross-linking does 
not impair the ability of dopamine to activate the D2AR, 
demonstrating that the receptor can bind dopamine and activate 
Gi with a disulfide bond between the C168 residue in each 
subunit of the dimer [49]. These observations suggest that the 


active receptor state remains within the receptor dimers 
themselves. 

The possible role for a dimeric species in B2AR signaling 
was discovered when it was found that receptor ligands affect 
the monomer:dimer equilibrium [48]. Isoproterenol increased 
the relative amount of dimer by 45% and also protected the 
dimer from the disruptive effect of the TM6 peptide [48]. These 
results suggest that isoproterenol stabilizes the B2AR dimer 
even in the absence of G proteins, which is compatible with our 
model. 

For any particular redox environment, there necessarily 
occurs an equilibrium between disulfide bonds and free 
sulfhydryls, which are then free to be influenced by the 
environmental pH. If two receptors dimerize by forming a 
disulfide bond, they will necessarily produce a pair of free 
sulfhydryls, because the four cysteine residues of the two 
disulfide bonds form only a single disulfide bond, thereby 
leaving two sulfhydryl groups free. This is a chemical 
disproportionation reaction that happens when a substance 
can be simultaneously oxidized and reduced, thereby allowing 
the formation of free sulfhydryl groups. 

We tested an extension of our previous model [16] in the larger 
context of past modelling efforts. Although this B2AR model 
without loops, in vacuum, with backbone restraints may not 
completely validate a GPCR activation hypothesis, our model 
offers a potentially powerful approach to modelling two-state 
interactions as the differences in interaction energies. These 
simulation results provide new insights into the relationships 
observed between selected biochemical data and the structural 
interactions that occur during receptor activation. 

In summary, these findings place the activation of the BZAR 
into a dynamic and general molecular model that explains many 
aspects of how ligands or mutations alter the B2AR to produce 
activation. These findings generally support Kenakin’s concept 
of protean agonism through a mechanism that describes ligand 
biased conformational change and suggest that there are multiple 
ways to activate and allosterically modulate the B2AR [2,50]. 
This model also relates the transmembrane domain motions of 
GPCRs to explicit receptor states and demonstrates that many 
elements from several different models may be important parts of 
a bigger view of receptor activation. Hopefully, this work will 
cultivate an appreciation for receptor function in terms of a 
general, two-state acid—base model that will present fruitful 
alternatives for the modeling of GPCR activation in the future. 


5. Conclusions 
The major conclusions from this work are: 


(1) The electrostatic interaction energy changes correlate well 
with the experimentally observed ligand efficacies. 

(2) The interaction energy changes for the helical/transmem- 
brane domain interactions are generally compatible with 
previously described events for GPCR activation. 

(3) This molecular model ties together many experimental and 
theoretical findings observed for the B2AR and other 
GPCRs. 
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